Noise Effects on Synchronized Globally Coupled Oscillators 
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The synchronized phase of globally coupled nonlinear oscillators subject to noise fluctuations 
is studied by means of a new analytical approach able to tackle general couplings, nonlinearities, 
and noise temporal correlations. Our results show that the interplay between coupling and noise 
modifies the effective frequency of the system in a non trivial way. Whereas for linear couplings the 
effect of noise is always to increase the effective frequency, for nonlinear couplings the noise influence 
is shown to be positive or negative depending on the problem parameters. Possible experimental 
verification of the results is discussed. 



Systems of coupled nonlinear oscillators are a generic paradigm of a whole class of problems arising in physics, 
chemistry, or biology Examples are Josephson Junction arrays (JJA) charge-density waves ||, thin film 
fabrication Q, chemical reactions and cardiac tissue ||, neuronal activity 0, and many more (see, e.g., Ref. B and 
references therein) . Due to the diversity of individual oscillators or to external (thermal) noise, these systems generally 
exhibit cooperative dynamical response or incoherent behavior as a function of the relevant parameters. The study 
of the transition between both regimes and the nature of the so-called synchronized phase requires the calculation 
of the phase distribution probability density. This involves solving a highly nonlinear, partial differential (Liouville) 
equation which is almost always a formidable task and only perturbative results can be obtained for practically all 
models of interest. 



In this paper, we approach these problems by introducing a new method, from a completely different viewpoint, 
that avoids dealing with partial differential equations at all. We focus on the stochastic (Langevin) evolution equation 
£j ' that allows to perturbatively obtain accurate results in a very simple and direct manner. This method can be applied 
to any system of N coupled homogeneous oscillators evolving in the presence of random forces, i.e., systems of the 
form 

J> 1 N 

<P t =0J- /(&) + ^ E - + <"*(*)> W 

where is the phase of oscillator i, lj > is the oscillator frequency, n is a constant, /(</>) is any 2-7r-periodic function, 
00 ■ T is any separable (see below) function, and the random term ?y, is local (uncorrelated from site to site), Gaussian 
^\ , and stationary. Aside from these requirements, r\i can be any process defined by a stochastic differential equation, 
e.g., it can be white or colored noise. To illustrate our method, we consider an Ornstein-Uhlenbeck process, given by 
i)i = -7Vi + 7&(*)i witn being uncorrelated Gaussian white noises [(&(*)&(£')) = &i 3 8{t - t')\ (• • •) stands for 
averages over noise realizations] and r — I/7 being the correlation time. 

Far from being academic, model (Q) contains already several important applications. To begin with, the coupling 
can describe many different systems: For instance, in the case of the phase approximation of globally coupled nonlinear 
oscillators or of a JJA, F(0j — 4>j) — sm (<& ~~ 4>j + a )i a being a constant, whereas in reaction-diffusion (RD) problems 
£j [ such as those arising in growth models ||4j , the coupling comes from the discretization of some spatial linear operator 
■ followed by a mean field approximation, yielding T{4>i — 4>j) = 4>i ~ 4>j- Second, the nonlinear term, /(</>), for the 
oscillator can be chosen as required in each problem; as in most coupled oscillator models, we take here f(4>) = sin</>, 
. but our method can be applied to other choices as well. 

We now begin the study of Eq. ([!]) in the case when, as mentioned above, T is separable, i. e., T{4>j — (pi) — 
^2 A n (4>i)B n (4>j). This technicality is needed to simplify the study but does not not restrict much the applicability 
of our calculations. Eq. (|l|) becomes then 



= u,-/(&) + -^5>"(&) 



A' 



ar,i(t). (2) 



The first step is to take the limit N — > 00. In this limit the sum within the brackets above can be computed in 
terms of the mean value of (f> As way of example, we focus on the RD problems and JJA's, for which we arrive 
(respectively) at 
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4> — lo — asin0 + k [((f>) — <f>] + 077(f), (3) 
cf> = lu — a sin <p + 

+k [cos(4> — a) (sin <f>) — sin(</> — a) (cos (f>)] + <J7](t) . (4) 

We note that we have thus reduced the N equations in (jl]) to a single (self-consistent) one. 

We deal first with eq. (^) for RD problems, hereafter called linear coupling model. When a = 0, it can be 
straightforwardly shown that if to < loc = a there is a stable time independent pinned solution <j) = sin~ cv/a, 
whereas if to > u>c the phase increases oscillatorily in time with frequency Q 2 = ui — a 2 . The natural question to ask 
is whether this scenario, involving a depinning transition at up, changes when noise is switched on, and if so, how. 
To address this issue we choose (0) as our order parameter, because in the pinned phase (0) is constant whereas in 
the depinned phase it increases in time. For noise intensities which are small compared to u) or k, we expand cue and 
(f> in powers of a, u c = a + a^a + a^a 2 + C(cr 3 ) and <p = tj>^ + ct0 (1) + a 2 ^ + 0(a 3 ). We first look for possible 
changes in wg, by inserting these two expressions in Eq. (|^); collecting powers of a and imposing that (0«) = for 
i = 0, 1, . . . (i.e., that the oscillators remained pinned), we can compute the a^' corrections to use- The final result is 

u>c -, o- 2 



a 4k(kt + 1) 



32k 4 [(kt+1) 



O(t) 



0(a 6 ). (5) 



From this expression, we immediately see first, that the transition occurs at a value of to which is lower than in the 
deterministic case, which means that the noise is actually helping the oscillators to overcome the barrier and start 
their motion. Another important conclusion is that for a given, fixed set of the other parameters, r controls the state 
of the system: whereas for low r values, the oscillators are depinned, as r increases the pinned phase is set on. We 
thus see the importance of the correlation time in the critical region of the system. 

Having found the changes in wc, we now turn to the dynamics in the depinned phase, when to > u>c- We again 
expand <f> in powers of cr and write down equations for each contribution, which read 

(<i>W) = u - asm{4>W); (6a) 

(0«) = 0; (6b) 

((0«) 2 ) = -2((0«) 2 )(k + acos(0(°))) + 2C(t), (6c) 

C = |-C(« + 7 + «cos(0W)); (6d) 

(0( 2 >) = -a(4> {2) ) cos<0<°>) + a (( ^ ))2) Bin<^°>). (6e) 

where ( = {(j^rf). Eq. ( |6e] ) implies that (tfi^) grows in time faster than (<j)^), and therefore our expansion is not 
correct. This problem, very well known in deterministic equations P|, can be cured by realizing that the system 
dynamics does not depend only on one time scale but on two, t and t* = fi(a)t. Among the different approaches one 
can use to deal with this, we choose Linstedt's method: we assume that <fffl depend on time through the combination 
t + t*; as for ju(cr), dimensional analysis shows that the noise term is of order one when ju(cr) = vu 2 . Introducing 
this new time scale in equations ( |6a[]6e[ ) and imposing the usual solvability condition on equation ( |6e[ ) (Fredholm's 
alternative, see e. g. for details) we obtain 



a I ' ((jyj) sin(0 o ) 



■dt, (7) 



2T 

where T = 2n/fl is the period of (cp -*) computed from eq. (6a). With this value, the effective frequency of the 



oscillations, defined as u e g — limt )S _ >00 ^[((f>(t + s)) — (<p(t))], is: 

lu cS = n(l + isa 2 +C(cr 4 )) (8) 
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FIG. 1. Predictions given by eqns. (ph up to order a 2 (solid line) and (^|) 
numerical solution of eqns. (^) for N = 5 000, a = 1, K = 1, u) = 1.1 and r 
calculation was done with an explicit 2.0 order weak scheme ju^]. 



up to order a 2 (dashed line) are shown with the 
= (squares) or r = 10 (circles). The numerical 



Figure [l] compares our analytical, order a 2 expression for Ct; e ff with the numerical solution of Eqs. (||) up to 
a ~ 0.5, already a not so small value. The expression for lo c q above indicates that the noise increases the frequency as 
compared with the deterministic value, in agreement with our previous conclusion that the noise helps the oscillators 
jump over the potential barrier: The larger the noise strength, the easier the potential barrier to overcome, effectively 
suppresing (renormalizing) it in the a — > oo limit. The fact that v > 0, implying that the frequency increases with 
the noise strength, can be proven rigorously when r = 0, i. e., for white noise, but we conjecture that this result is 
general for linearly coupled systems. 

Another advantage of our approach is that it can be applied to other limits, such as the case when the noise intensity 
is larger than u> and a. Expanding the solution now in powers of a, the same steps as before yield, up to order a 2 , 



a 2 



exp 



a 2 (KTe 



— s/ T 



2k(kt - 1) 



K 2 L0 2 + UJ^ 

ks > sinws ds 



(9) 



where a = (j/^/kt 
very accurate, in this case for values of a down to a 



1. Again, the comparison with the numerical solution (see Fig. shows that our approach is 

1.2, very close to the values of ui and a. We want to stress that, 
except for a small interval 0.5 < a < 1.2, our analytical predictions (@) and (||) describe to a high degree of accuracy 
the numerical solution for any noise intensity, providing a global picture of the main features of the linearly-coupled 
system behavior. 
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FIG. 2. Comparison of the numerical values of w e ff for the linear (squares) and KS (circles) models with the small a 
approximation (g) (solid lines). Parameters are as in Fig. [ljwith r = 0. Note the synchronization-desynchronization transition 
for the KS model is at a ~ 0.5. The inset shows the computed time evolution of the order parameter, p, (see text) in the KS 
model for a = 0.2 (circles) and the comparison with our analytical aproximation (solid line). 



We now move on to the nonlinear problem (|J) in the case a = 0; Eqs. ([j]) are then known as Kuramoto-Sakaguchi 
model (KS). In the small noise regime, our technique applied to the KS model leads once again to Eqs. ( |fo|]6e| ). This 
stems from the fact that, when a — > 0, the phases of individual oscillators are similar, and consequently sm((/)j — 
fa) ps (f>j ~ fa. However, due to the periodicity of the KS coupling term, the model is invariant under changes 
fa — (f)j <-> fa — (f>j + 2mr, for any integer n. Therefore, 2tt jumps between oscillators are possible without consequence 
other than the breakdown of the linear approximation. Of course, the larger the noise strength, the more likely such 
jumps are, and the KS model becomes desynchronized at a finite value of a, i. e., it has a true synchronization- 
desynchronization transition at a 2 = n when a = jlj, and at lower values for a ^ Up to that point, it can 
be seen from Fig. || that the linear model, the KS model and our analytical prediction Eq. (|^) are all in excellent 
agreement with each other. For our parameters, a = 1, u> = 1.1, the transition takes place at a sa 0.5. 




FIG. 3. Diagram of the possible signs for v(a) for different values. Leftmost line corresponds to r = 10, rightmost one to 
t = 0. In both cases, Region I (y a > 0) is ond the left of the line, Region II [y a < 0) lies on the right. The square and the 
circle points correspond to the cases analyzed in figure (U|). 
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The order parameter for the synchronization-desynchronization transition is p = (cos < 
chronized phase, p oscillates around a nonzero value, while in the desynchronized phase, p 



) 2 ; in the syn- 
0. The inset in Fig. | 



compares the numerical value of r with the result of our aproximation, p = 1 — a 2 ((</>W) 2 ) + 0(a 3 ), exhibiting the 
quantitative validity of our approach. The agreement is very good up to the critical value of a. As for the limit 
a — ► oo, the computations, although feasible, are quite involved, and therefore we do not present the results here. 

In the two cases analyzed so far, in the synchronized phase noise helps the oscillators overcome the nonlinear 
potential, hence effectively increasing their frequency. However, this intuitively reasonable picture is not the generic 
situation: Antisymmetric couplings verifying T(—<p) — —T(<f>) keep the difference between oscillators small, forcing 
their motion to be approximately the same, and then the above picture holds, but if the coupling lacks this symmetry 
(like V(4>i — 4>j +a) for nonzero a), the oscillators do not tend to synchronize (they rather have phases separated by a 
factor a). In the latter case, the noise competes with the frustration induced by the coupling and its activation effect 
disappears. Our calculations for Eqs. (Q) when a ^ and a — > show that this is indeed what occurs. Our method 
leads now to a value of u, the a 2 correction to the frequency, which depends on a, given by 



1 

2T 



((f) 2 ) (a sin(0 o ) — 2k sin a) 
luq — asin(0o) 



dt 



(10) 



where usq = us + Ksina, this parameter playing the role of an effective driving. Eq. ( pL0| ) shows that, when a < ao = 
sin~ (1/2k), v a is positive or negative depending on the value of us, whereas when a > ao we have v a < 0. This 
divides the (u,a)-space in two regions, according to the possible signs of u a , as depicted in Fig. ||[ Interestingly, in 
region I a (or equivalently us) acts as a switching parameter between noise induced acceleration or deceleration of 
the oscillator motion, as illustrated in Fig. ^. Note that this is a novel, noise intrinsic phenomenon, absent in the 
deterministic problem, whose origin is the nonlinear coupling: If we add a constant to the coupling in Eq. (||) only 
the external frequency is changed (us — > us + na) . 
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FIG. 4. us c s obtained from simulations for uio = 1.006, r = 0, ft = 1 and a = 0.5 (squares) or a = 0.55 (circles), compared to 
the analytical approximations (solid lines) given by eqns. (^) and (|To[) . [Squares and circles refer to the corresponding points 
in Fig. |. 



In summary, we have introduced a procedure to study systems of globally coupled oscillators in the synchronized 
phase, general enough to analyze any stationary random processes generated by a stochastic differential equation, 
any form of the nonlincarity, and a large family of coupling terms. We have used it to analyze problems relevant in a 
number of physical contexts, obtaining fairly good results in a much more straightforward way than the approaches 
proposed so far. For linear coupling, we have computed the noise-induced corrections to the critical value for the 
depinning transition as well as the effective frequency in the depinned phase for almost all values of the noise. We 
have shown that the same results hold for the KS model. In the general case of frustrated JJA's we have found 
that the effect of noise on the effective frequency is not trivial, leading to an increase or a decrease of the effective 
frequency depending on the system parameters and the driving frequency us. This very important result shows the 
power of our technique as it was not known prior to our work. Our predictions can be directly checked by experiments 
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in very many fields, among which we suggest JJA's subject to a magnetic field where tuning the voltage for a given 
temperature one should be able to find the two different behaviors predicted. 

We are indebted to J. A. Cuesta and R. Cuerno for discussions. Work at GISC (Leganes) was supported by CICyT 
(Spain) grant MAT95-0325 and DGES (Spain) grant PB96-0119. 
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